Quantum Monte Carlo study of a one-dimensional phase-fluctuating condensate in a 

harmonic trap 
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We study numerically the low-temperature behavior of a one-dimensional Bose gas trapped in 
an optical lattice. For a sufficient number of particles and weak repulsive interactions, we find a 
clear regime of temperatures where density fluctuations are negligible but phase fluctuations are 
considerable, i.e. , a quasicondensate. In the weakly interacting limit, our results are in very good 
agreement with those obtained using a mean-field approximation. In coupling regimes beyond 
the validity of mean-field approaches, a phase-fluctuating condensate also appears, but the phase- 
correlation properties are qualitatively different. It is shown that quantum depletion plays an 
important role. 

PACS numbers: 03.75.Hh, 03.75.Lm,05.30.Jp 



In a spatially homogeneous one-dimensional (ID) gas 
of bosons, spontaneous symmetry breaking is excluded at 
all temperatures T P, 0, Q . Long-range order is absent 
in this system due to the fluctuations of the phase of the 
order parameter, as seen in the asymptotic decay of the 
equal-time single-particle Green's function (algebraic for 
T = and exponential for T > 0) 0, H H 0- However > 
analytical studies of trapped ID bose gases with a fixed 
particle number, N, reveal new phenomena. For exam- 
ple, the ground state of a nonintcracting ID bose gas in 
a harmonic trap with frequency u> becomes macroscopi- 
cally populated below a temperature T c ~ Nhui/ ln(27V), 
i.e. , a Bose-Einstein condensate (BEC) exists in the fi- 
nite system Q. Several mean- field studies indicate that 
weak repulsive interactions introduce an additional effect, 
namely, that fluctuations of the density are suppressed at 
low temperatures, and a phase-fluctuatin g co ndensate, or 
quasicondensate, appears & 1^, 111. 113, Tl3l. 14 1; for the 
homogeneous system see [15j. Roughly speaking, this is 
a regime in temperature where the distribution of par- 
ticles in space is given by a temperature-independent 
Thomas-Fermi profile, while thermal fluctuations of the 
phase are present and lead to a phase-coherence length 
that is smaller than the condensate cloud. Only below 
a much lower temperature do phase fluctuations also be- 
come negligible and phase coherence extends over the 
complete condensate cloud. A phase-fluctuating conden- 
sate emerges only in low dimensions and has been experi- 
mentally observed in sufficiently anisotropic 3D trapping 
geometries flij . 

In this paper, we verify the existence of a ID quasi- 
condensate, on trapped optical lattices, starting from a 
microscopic approach whose validity is not restricted to 
certain parameter regimes. Using quantum-Monte Carlo 
simulations of the Bosc-Hubbard model, we investigate 
the properties of the phase-fluctuating condensate for 
various choices of interaction strength and density. For 



weak interactions, our results are in excellent agreement 
with mean-field estimates in the continuum [lOfl . In an 
intermediate-coupling regime beyond the weakly inter- 
acting limit, but not yet in the strong-coupling domain, 
the quasicondensate regime spreads over an even larger 
temperature range. Nonetheless, we observe qualitatively 
different phase-correlation properties which do not follow 
from existing analytical approaches. 

First, we briefly review the notion of a quasiconden- 
sate (QC) in a weakly interacting (WI) ID trapped bose 
gas as presented in [l(J. In the WI limit, we have 
7 = mg/h 2 n <C 1, where n is the average particle den- 
sity, m the particle mass, and g the coupling constant of 
a repulsive contact interaction 17|. In second-quantized 
representation, the Hamiltonian of the system is given by 



H = J dzft(z) 



(1) 

where ip(z) is the bosonic field operator and V t (z) = 
muj 2 z 2 /2 the external trapping potential centered at the 
origin. The authors in [lfj show that density fluctua- 
tions (Sn(z)Sh(z')) , where n(z) = n(z) + 8n(z), are sup- 
pressed for inverse temperatures (3 3> Pd = {Nftio) -1 
(Boltzmann constant ks = 1). The field operator can 
then be expressed as ip(z) = y n(z) cxp[i0(z)] . Thermal 
fluctuations of the phase are evident from the decay of 
the one-particle density matrix, or (equal-time) Green's 
function, which is obtained as [lfj 



(^(z)iP(z')) = vX^Mz')cxp(-(^L')/2), 
where 54> 2 ZZ , = {<j>{z) — 0(z')) 2 , and 
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where /xtf = (3iV g / '4) 2 ^ 3 (muj 2 /2) 1 / 3 is the chemical po- 
tential, and Ltf = V^Mtf jmuj 2 the half size of the 
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FIG. 1: (Color online)(i) Thermal regime: As a result of both 
density and phase fluctuations, the widths of the density pro- 
file, L n , and the momentum profile, L k , depend on j3 with 
L n (f3) > L^ at , L k (P) > LJ at . (ii) Quasicondensate: Den- 
sity fluctuations are negligible (L„ = L^ at , L is the half size 
of the condensate cloud for /3 > /3qc). However, the phase 
fluctuates, as seen for Lk(f3) > L% at , phase correlation func- 
tion C(0,0.8L) < 1. (iii) True condensate: Phase fluctua- 
tions also disappear as Lk(P) approaches -LJ at , and C(0, 0.8L) 
approaches 1. The numerical results are /3qc ~ 0.03 and 
/3tc = 0.115(5) (analytical estimate: /3tc = 0.116). 



condensate in the Thomas-Fermi (TF) approximation 
(mtf ^ see e.g. [3]). For (3 — ► oo, it follows from Eq. 
(pE|) that (S(j> 2 zz ,) -> 0, and thus G{z,z') -> v^Mz'), 
i.e. the system is completely phase coherent. 

As is well known [18[, the discretization of Eq. ((T|) 
yields the Bose-Hubbard model with lattice Hamiltonian 

h= -jj2 (fij-fij' + H - c ) + 1 E «j - 1) 



-V t £(j-M/2) 2 



(3) 



where creates a particle at optical lattice site j and 

hj = OjCij. The number of lattice sites M is chosen 
such that the occupation at the boundaries of the lat- 
tice is zero. We work in units where the lattice spacing 
and H 2 /2m are set to 1. In these units, parameters in 
Eqs. (p} and (J3j) are related as follows: J = 1, U = g, 
Vt = (fno) 2 /4. Our system is defined on an optical lattice, 
however, in the WI and degenerate limits, the discrete 
and continuum description of the phase-fluctuating con- 
densate arc equivalent 



12] 



We investigate this model 
using a worm update quantum Monte Carlo (QMC) 
method in the canonical ensemble [20| . This nonlocal up- 
date scheme allows for an efficient evaluation of the equal- 
time Green's function G(j,f) = (aja,-'). Phase correla- 
tion properties are also apparent from the shapes of the 
momentum profile, n(k) = J2j j' G(j,j') exp[ik(j — j')], 
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FIG. 2: (Color online) N = 2000, U = 0.4, Vt = 0.3 (7 = 
0.003). Density fluctuations are suppressed for j3 > /3qc ^ 
0.03 and the density profiles iij = (hj) (o) are identical for 
P > /3qc (b)-(d). In (a), the line is a fit to a Gaussian profile, 
while in (b)-(d) the lines are fits to a TF profile no(l — j 2 /L 2 ) 
where L = 12.6 (Ltp = 12.6). Error bars in the figures are 
always smaller than the symbol size. 



and the rescaled Green's function (phase correlation func- 
tion), C(j,f) = G{j,j')/y/njfij7, where nj = (hj) is the 
density distribution. In the WI, TF and mean-field (N 
large enough) limits, we observe three different regimes 
in temperature: the thermal regime (TR), the quasicon- 
densate (QC) and the "true condensate" (TC). These 
regimes are separated by smooth crossovers, which are 
characterized by the temperatures 1//3qc and 1//?tc- 
We discuss the properties of the three regimes and com- 
pare our results with the mean-field results listed above. 
Furthermore, we consider the emergence of the phase- 
fluctuating condensate depending on the choice of pa- 
rameters N, U and Vt, including parameter sets that are 
beyond the WI limit. 

In the thermal regime (3 < /3qc, both the density and 
the phase are governed by thermal fluctuations. Hence, 
all quantities exhibit a strong temperature dependancc. 
The width of the density profile, L n (f3) (standard devia- 
tion of rij), decreases throughout the TR, until it reaches 
its minimum value L^ at at the inverse temperature (3qc, 
as shown in Fig. [TJ The shape of the density profile is ap- 
proximately Gaussian, which is expected for high temper- 
atures where the bosonic nature of the particles becomes 
less relevant [Fig. [5] (a)]. The strong phase fluctuations 
are seen in the width of the momentum profile, Lk(f3) 
[standard deviation of n(k)], which is much larger than 
its minimum value L| at [Fig. [1] , as well as in the exponen- 
tial decay of the Green's function [Fig. [3] (a)]. We recall 
that there exists no analytical estimate for /3qc; density 
fluctutations are merely predicted to become small for 
[3 > f3 d = (AVtN 2 )- 1 / 2 [lj|. For the parameter set in 
the Figs. QI3 , we have j3 d w 0.02/3 Qc . 

In the quasicondensate (/3qc < P < Ptc), th e density 
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FIG. 3: (Color online) N = 2000, U = 0.4, Vt = 0.3. 
The density profiles *Jnyn~j (o) are virtually independent 
of temperature for j3 > /3qc ~ 0.03 (b)-(d), while the 
Green's function G(0,j) (A) varies with (3. The lines 
are the analytical estimates, a Gaussian fit for y/nonj in 

(a), a TF profile no\/l — j 2 /L 2 for ^/norij in (b)-(d) and 
^n5nJe X p(-a|ln[(L-i)/(L + i)]|) [see Eq. (2)] for G(0,j). 



no longer fluctuates: the density profiles nj at differ- 
ent temperatures [3 > /3qc are identical and in excellent 
agreement with a TF inverse parabola shape of half size 
L, where L equals to Ltf [Figs. EJb)- E^d)]. However, 
thermal fluctuations of the phase are considerable. Thus, 
the width of the momentum profile, and the phase corre- 
lation function still change with decreasing temperature 
(Fig.[T]). Figures[3jb)-[3]Jd) demonstrate that the Green's 
function varies substantially for different j3 > /3qc , while 
the density profile ^n^rij is invariant. We find very good 
agreement with the mean- field result Eq. (2): the ansatz 
^/riQUj cxp{— a\ ln[(L — j)/(L with a a fitting pa- 

rameter, reproduces G(0,j) (Fig. [3]). 

At temperatures much lower than 1//3qc, phase fluc- 
tuations also become suppressed. Therefore, it is 
meaningful to introduce a second crossover tempera- 
ture, 1//?tC) to a phase-coherent regime, the true con- 
densate. We define /3tc as the temperature where 
C(0,0.8L) = exp(— 0.25). This definition is somewhat 
arbitrary, but yields a scale below which the system 
can safely be considered to be phase-coherent. In ad- 
dition, it allows for a comparison of analytical and nu- 
merical results, with the analytical equivalent of /3tc 
being (3tc = ^ftd^TF ln(9)/37kj [which is the temper- 
ature where {5(j^ z ,) = 0.5 for \z'\ = 0.8Ltf; see Eq. 
2]. Note that the true condensate is not to be con- 
fused with a BEC, and strictly only appears at zero 
temperature (Thomas- Fermi condensate; see [13]). In 
the true condensate, the Green's function and the mo- 
mentum profile approach their temperature-independent 
ground state shapes, as shown by the convergence of 
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FIG. 4: (Color online) The size s = /3tc//3qc, of the QC 
regime depending on the number of particles N (a - c, U = 
0.5, V t = 0.01) and coupling constant U (d - f, N = 2000, 
V t = 0.4). 

L n /LT (o), L k /LT (A), C(0,0.8L) ( V ) 

a) 7 = 0.08, /3 TC = 3.5, /3 TC « 3.2(4), /3 QC « 2.0, s « 1.6 

b) 7 = 0.04, p TC = 2.5, (3 TC « 2.6(4), f3 QC « 1.0, s « 2.6 

c) 7 = 0.02, p TC = 1.64, p TC « 1.6(1), (3 QC ~ 0.3 , s « 5.3 

d) 7 = 0.001, /3 TC = 0.06, (3 TC » 0.06(2), /? QC w 0.04, s « 1.5 

e) 7 = 0.006, Arc = 0.15, Arc « 0.15(5), /3 QC ~ 0.03, s « 5 

f) 7 = 0.17, /3tc = 0.82, NO fe, Pqc « 0.01, sw^ = 100 
Comparison of (a)-(c), and (d)-(f), respectively, shows that s 
increases with increasing N and f7. 



L k ((3) and C(0,0.8L) in Fig. [T] The system becomes 
practically phase coherencent, as illustrated in Fig. (3Jd), 
where ^nonj m G(0, j). For the parameter set in Figs.[I|- 
121 we find that /3tc ~ 0.115(5) which agrees with the an- 
alytical estimate /3 T c = [{U 2 /6NV 2 )^ 3 ln(9)] = 0.116. 

We now study the appearance of a phase-fluctuating 
condensate depending on the system parameters U, Vt, 
and N. The WI limit is characterized by 7 = U/2n <C 1. 
We define the average density by n — N/2L. Since we 
find that L equals Ltf within error bars in all cases, we 
use 7 = (WyWtN 2 ) 1 / 3 . Note that the magnitude of 7 
in the inhomogenous system differs from that in the ho- 
mogeneous system. In Figs. QJa)- 0|c), we demonstrate 
the effect of varying the number of particles N at on-site 
repulsion U = 0.5 and trapping Vt = 0.01. Both (3qc and 
/3tc decrease with increasing N, but since the change of 
Pqc is much greater than that of /3tCj the region of the 
phase-fluctuating condensate, s = /3tc//3qc, increases 
with increasing N. The effect of varying U is illustrated 
in Figs. Btd)-Stf), where N = 2000 and V t = 0.4. In- 
creasing U, causes (3qc to decrease, while /3tc increases. 
Thus /3tc//?qc grows with increasing U. More generally, 
the deeper we are in the Thomas-Fermi limit of large 
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FIG. 5: (Color online) N = 2000, U = 10.0, Vt = 0.4 
(7 = 0.17). The density profile y/JTonj (o) is invariant for 
P > /9qc ~ 0.1, while the Green's function G(0,j) varies 
throughout the QC regime. However, G(0,j) does not ap- 
proach y/noWj (as in Fig. The line is a TF fit to ^/nonj 
with L = 33.3 (Ltf = 33.5). 




j j 

FIG. 6: (Color online) Ground state y / nonJ (o) and G(0, j) 
(A) for different 7. The lines are fits to a TF profile with (a) 
L = 35.2 (Ltf = 35.6) and (b) L = 48.0 (Ltf = 48.3). For 
increasing 7, the correlation hole becomes more pronounced, 
and thus the quantum depletion increases. 

fi TF /huj ~ [(NU) A /V t ] 1,& , the larger is the size of the 
quasicondensate. 

While for the parameter sets in Figs. 01a)- 0|e), good 
agreement of G(0, j) with expression (2) is observed for 
all /3, as well as /3tc ~ Ptc, this does not apply to the 
example in Fig.0Jf), where 7 = 0.17. Instead, G(0,j) ex- 
hibits a qualitatively different behavior when approach- 
ing the ground state, as can be seen in Fig. [5] the phasc- 
nuctutating condensate is beautifully realized, however, 
G(0,j) does not approach ^Jn^rij for T — ► 0, and cannot 
be described by Eq. (2). 

We consider the effect of stronger interactions in the 
trapped system on the phase-correlation properties in 
more detail. In Fig. [6l we show the ground state pro- 
files for N = 300, V t = 0.01, U = 2.0 and U = 5.0, 
respectively. We observe the same qualitative behaviour 
of the saturated Green's function as in Fig. Close to 
the center of the cloud, where the density is higher, the 
decay of G(0, j) is exponential; however, it broadens to- 



ward the outer, more diluted regions of the cloud. By 
comparing Figs. [6ja) and [6jb), it can be seen that for 
larger 7, the correlation hole in the central region of the 
condensate is larger. Clearly, the quantum depletion for 
the parameter sets in Fig.[B]is much more significant than 
for systems in the WI limit, since the number of particles 
with zero momentum N(k = 0) = J2j j' G(j,j')/M de- 
creases if G(j,f) < y/nJnJT. We also observe this behav- 
ior if the density is smaller than 1 everywhere in the trap, 
and therefore exclude the possibility that it is an effect 
of the optical lattice. The shape of G(0, j) is consistent 
with the two limiting cases, i.e. exponential decay in the 
strong-coupling limit, and phase-coherence [i.e. Eq.(2)] 
in the weak-coupling limit. 

We conclude with a summary of our main results. In 
the mean-field, weakly interacting and Thomas-Fermi 
limits, both true a Thomas-Fermi condensate and a 
phase-fluctuating condensate emerge. Phase correla- 
tion properties, manifest in the characteristic decay of 
the single-particle density matrix, agree with surpris- 
ing precision with the mean- field theory in [lo| . In 
an intermediate-coupling regime, a true condensate no 
longer appears. The regime of the phase- fluctuating con- 
densate persists even longer, extending to T = 0. We ob- 
serve a qualitatively different decay of the Green's func- 
tion, which cannot be accounted for by mean-field stud- 
ies. 
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